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SUMMARY 


A procedure is described for improving a linear nonproportionally damped 
analytical model of a structure. An approximate analytical model is assumed 
to be available. Given also is an incomplete set of measured frequencies, 
damping ratios and complex mode shapes of the structure, as may be obtained 
from a vibration test. A method is developed which finds the smallest changes 
in the analytical model such that the improved model matches the measured 
modal parameters. This report includes a discussion of the problem, the deri- 
vation of the algorithms and a description of the computer implementation. An 
example of an application of the procedure to improve an analytical model is 
presented. 
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INTRODUCTION 


With the increasing complexity and accuracy requirements of modern aero- 
space structures, system identification may be perceived as an attractive 
alternative to intuitive validation of analytical models. In recent years, 
within the framework of linear undamped systems, many procedures have been 
proposed to improve an analytical model by using measured modal parameters. 
Typical publications which are related to the method developed in this report 
are represented by Refs. 1-7. 

Analytical mass and stiffness matrices are often available from some 
source, such as a finite element analysis. Certainly, the resulting model is 
only an approximate representation to the structure and contains errors of 
undetermined magnitude. On the other hand, vibration tests result in incom- 
plete sets of natural frequencies and mode shapes of the structure. These 
data also contain errors whose magnitudes are unknown. Assuming it is desired 
to have an analytical model which will be consistent with the analysis and the 

test, Berman and Wei^ applied the Lagrange multiplier method of Baruch and Bar 
3 

Itzhack to minimize a weighted matrix norm and obtained an improved model 
which exactly predicted the measured modal data and also minimized the changes 
in the analytical model. 

The acceptability of the original analytical model as a baseline for fur- 
ther analysis depends upon the size of the minimum changes. The requirement 
for a small modification implies that the analytical model is a good represen- 
tation of the structure. Evaluation of these changes is based on the user's 
judgement and statistical information regarding the changes. In contrast with 
more conventional procedures which require trial-and-error adjustment of the 
oriqinal model to match the test data, this method yields an improved model 
and information of possible deficiencies in the test data or analytical model. 
Applications of this procedure to improve the analytical models of large 
structural systems have been reported**^. 
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The procedure mentioned above, among others, is based on the assumption 
that the system possesses only negligible damping or that the damping is pro- 
portional and that the measured modal parameters are real. However, there are 
cases in which damping effects must be taken into account. For instance, 
ground vibration tests of a helicopter fuselage have shown that the modal 
damping ratios of several lowest modes are over 12 percent and the mode shapes 

O 

are essentially complex . Therefore, it was necessary to refine the previous 
method so that the more realistic, nonproportional ly damped system and the 
corresponding modal parameters may be equally treated. In the analyses repor- 
ted hereafter, the system damping can be of any type so long as it can be rep- 
resented by a real symmetric matrix coefficient of the velocity vector in the 
dynamic equation. 

Primary concerns in the present study are as follows: 

1. The method must be applicable to large structural systems. Let the 
number of degrees of freedom be n and the number of measured modes be m. Typ- 
ically, n is of the order of hundreds or thousands and m is much less than n. 
Methods involving extensive iteration or repeated eigensolutions may be prohi- 
bitively expensive and the convergence problem would always involve some 

2 

uncertainty. Methods requiring the solution of sets of equations of order n 
would also be impractical. 

2. Besides the analytical mass (M^) and stiffness (K^) matrices, the 
present procedure requires an analytical damping matrix (C ^). Each of these 
analytical matrices is to be minimally changed to agree with the test data. 
It is acknowledged that no general rules currently exist for formulating a 
damping matrix as are those for mass and stiffness matrices. One way to 
obtain in the first approximation is to represent it as a proportional 
matrix which is a linear combination of M^ and K^. This step, of course, 
involves considerable engineering judgement and the accuracy may not be high. 
It should be noted that the procedure will modify this matrix to make it 
appropriately nonproportional if the modes are complex. 
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3. Unlike its undamped counterpart, a damped system has complex-valued 
modal matrices n and <j>. In order to guarantee that the improved matrices are 
real symmetric, the constraints which must be imposed on the system equations 
render the problem considerably more complicated than would otherwise be the 
case. 


A method was developed which is applicable to realistically large models, 
is computationally efficient, and satisfies all the theoretical requirements, 
except that further improvement of the identified stiffness matrix is 
required. Further research is necessary to eliminate this deficiency. The 
modified analytical model is, however, a distinct improvement over the orig- 
inal . 
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SYMBOLS 


A 

B 

C 

E 

F 

f 

f 

t 

G 

g 

H 

h 

I 

K 

k 

l 

M 

m 

n 

P 

Q 

R 

S 

U 

Y 
W 
X 

Y 
Z 
a 
8 

Y 
A 


Matrix defined in Eq. (4) 

Matrix defined in Eq. (5) 

System damping matrix 
Coefficient matrix, Eq. (14) 

GH, defined in Eq. (27c) 

Matrix defined in Eq. (15) 

Matrix defined in Eq. (22c) 

Matrix defined in Eq. (27a) 

Matrix defined in Eq. (22a) 

Matrix defined in Eq. (27b) 

Matrix defined in Eq. (22b) 

Identity matrix 
System stiffness matrix 

w k kw k 

Number of test points 
System mass matrix 
Number of measured modes 
Number of degrees of freedom 
Matrix defined in Eq. (27d) 

Matrix defined in Eq. (27e) 

Real part of $ 

Imaginary part of $ 

Real part of 4>n 
Imaginary part of 
Weighting matrix 
Displacement vector, Eq. (1) 

Vector of system state variables, Eq. (3) 

Matrix defined in Eq. (27f) 

Matrix consisting of Lagrange multipliers, Eq. (18) 

Matrix consisting of Lagrange multipliers, Eq. (18) 

Matrix consisting of Lagrange multipliers, Eq. (18) 

Matrix defined in Eq. (14b) 
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e Weighted matrix norm 

0 y - y t 

A Matrix consisting of Lagrange multipliers, 

$ Rectangular complex mode shape matrix 

♦ W" 1 * 

i|> Lagrangian function 

Diagonal complex matrix consisting 
of damping ratios and natural frequencies 


Superscripts 


(•) 

Differentiation with respect to time 

n 

Complex-conjugate of a matrix 

T 

Matrix transpose 

* 

Conjugate transpose of a matrix 

+ 

Moore-Penrose inverse of a matrix 

Subscripts 

A 

Analytical matrix 

C 

System damping 

I 

Imaginary part 

K 

System stiffness 

M 

System mass 

R 

Real part 


. (14b) 



THEORETICAL DEVELOPMENT 


Theoretical Background 


A linear structural system of n degrees of freedom is described by the 
matrix equation 

MX + CX + KX = 0 (1) 


where M, C and K are real symmetric constant coefficient matrices. With the 
state-space representation, Eq. (1) can be rewritten as 


AY + BY = 0 


where 



( 2 ) 

(3) 


A = 


0 

M 


M 

C 


( 4 ) 


B = 



0 

K 


(5) 


By the use of Eqs. (2) to (5), the following orthogonality relationships are 
obtained: 


rd>oi T r*oi 


[• ] a [# ] = 1 

(6) 

r*o"iT 


|TJ B [«] - - <1 

(7) 

or equivalently. 


*> T M $ n + n 4> t m * + C * = I 

(6a) 
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and 


ft* T M$ft-$ T K$ = ft (7a) 

where ft is a complex-valued diagonal matrix containing modal damping ratios 
(real part) and natural frequencies (imaginary part) and $ is the correspon- 
ding mode shape matrix. In this analysis $ is a rectangular matrix and a is 
diagonal, of order m. With these modal parameters, the dynamic equation (1) 
can be cast in the form of 

M $ fl 2 + C $ n + K t = 0 (8a) 

Given ft and $, it can be shown that if M, C and K satisfy Eqs. (8a) and 
(6) then they also satisfy Eq. (7). In other words, Eqs. (8a) and (6) imply 
Eq. (7), but Eqs. (8a) and (7) do not imply Eq. (6) unless ft is nonsingular. 
Therefore, Eqs. (8a) and (6) have been selected to improve the analytical 
model since it is desired to include rigid body modes if the system has free 
boundary conditions. 

Intuitively, it might seem appropriate to improve M, C and K by directly 
identifying A and B of Eqs. (4) and (5). Unfortunately, this cannot be done. 
The reason is that the optimization techniques would distribute variations to 
all variables involved and, for example, the upper-left portion of the modi- 
fied matrix will no longer to be precisely zero, which is certainly unaccep- 
table. Thus, it is necessary to take the more difficult analytical approach 
of specifically identifying the physical matrices. 

Full Mode Computation 

Since in practice, the degrees of freedom measured in a vibration test 
are only a small portion of those specified in the analytical model, an inter- 
polation procedure is employed to compute the modal displacements at the 
unmeasured degrees of freedom. This preprocessing procedure, referred to as 
"full-mode computation," has been discussed in Ref. 6. A slight modification 
is made here to include the system damping. 
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Let be the ith column of the mode shape matrix $. The portion of * 
whose elements are measured is denoted by and the rest by With this 
coordinate transformation, Eq. (8a) can be written as: 



where, for the ith mode, ^ is the measured eigenvalue and M, C and K are par- 
titioned as shown. At the start of the identification procedure, M, C and K 
are the analytical matrices M A , C A and K A . It is apparent from Eq. (8b) that 

*2i = " (ft i M 4 + Q i C 4 + K 4 )_1 (n i M 2 + fi i C 2 + K 2 } *11 (9) 

Because analytical matrices are involved in the equation, the accuracy of 
$2-,- solved from Eq. (9) is expected to compare to that of M^, C A , or K^, even 
if is exact. Note that even if ^ is inexact, the improved model should 
exactly predict , which is the only data actually measured. The first step 
in the procedure is to compute the full modes from Eq. (9). The following 
analysis uses the full, but incomplete, modal matrices $(nxm) and fi(mxm). 

Improvement of Mass and Damping Matrices 

The basic equation to be used for improving analytical mass and damping 
matrices is the orthogonality relationship, Eq. (6a). Since m < n, given * 
and n, there are infinite sets of M(nxn) and C(nxn) satisfying Eq. (6a). A 
unique set of mass and damping matrices will result if, in addition to satis- 
fying the orthogonality relationship, a norm 

e = ||H h (M-M a ) W M II + ||W C (C-C A ) W c ll (10) 

is also minimized, where W^(nxn) and Wj,(nxn) are arbitrary real symmetric, 
nonsingular weighting matrices. Let * = R + iS and = U + iV; then the real 
and imaginary parts of (6a) are separated into 
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(11a) 


r t mu + u t mr - s t mv - v t ms + r t cr - S T CS = I 


r t mv + v t mr + s t mu + u t ms + r t cs + s t cr = 0 


(lib) 


Define the Lagrangian function 


mm TTTTTT 

V = e + l l (A_). (R'MU + U 1 MR - S ' MV - V'mS + R 1 CR - S'CS - I).. 

i=l j=l R 1J 1J 


mm TTTTTT 

+ l l (A.).. (R MV + V MR + S MU + U MS + R CS + S CR) .. (12) 

i=l j=l 1 1J 1J 

where A R and A^ are Lagrangian multipliers. Differentiating f with respect to 
M and C, respectively, and equating the results to 0 yield 

M = M A - ! W ‘ 2 (RA r U T + UA r R T - SA R V T - VA R S T + RAjV 1 + VAjR 1 


+ SAjU 1 + UAjS T ) W ‘ 2 

c = C A - | W ' 2 (RA r R T - SA r S T + RAjS T + SAjR 1 ) w - 2 


(13a) 


(13b) 


Obviously, requiring A R and Aj to be symmetric in Eqs. (13a, b) is equi- 
valent to imposing the constraints of M and C to be symmetric in Eqs. (11a, 
lib). The Lagrange multipliers (A ).. and (A.)., may be solved numerically 

K lj 1 1 J 

from the following set of equations obtained by substituting Eqs. (13a, b) 
into Eqs. (11a, b). 


* E ll'lklj * E 12*iknj J*Vk»[ J ( Vij 

r'M 

J J v J 

which, for simplification, may be written as 


(14a) 


EA = A 


where 


A r = 2 (I - r t m a u - U T M ft R + s t m a v + v t m a s - r t c a r + s t c a s) 


(14b) 


(14c) 
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( 14d) 


Aj = 2 (r t m a v + v t m a r + s t m a u + u t m a s + r t c a s + s t c a r) 

-? -?. 

The submatrices E^, E^, E^j and ^22 are functions of R, S, U, V, and W^' 
and their explicit function forms are given in the Appendix. In Eq. (14), a 
repeated index in a product denotes summation over the range of the index. 

Since A R and Aj are symmetric, there are m(m + 1) unknowns along with 
m(m + 1) equations in Eq. (14). In the implementation of the technique, the 
coefficient matrix E which has m(m + 1) x m(m + 1) elements is formulated by a 
FORTRAN program. The improved mass and damping matrices, M and C, result when 
A r and Aj are solved from Eq. (14) and are numerically substituted into Eqs. 
(13a, b). 

In this section, a numerical method has been used to solve m(m + 1) equa- 
tions. In practice the number of measured modes is small enough to make this 
solution economically practical. 

Improvement of Stiffness Matrix 

Given the full modal matrices $ and n, the improved M and C, and the ana- 
lytical stiffness matrix K A> an improved stiffness matrix K must satisfy the 
dynamic equation 

<t> T k = - (0 2 $ T M + M T C) W k (15) 

or 

<|> T k = f (15a) 

and the constraints 

k = k (16a) 

k = k T (16b) 

in which an asterisk denotes the conjugate transpose of a matrix, <j> = Wj^ <t, 
k = W|^ K and (nxn) is an arbitrary, real symmetric, nonsingular weight- 
ing matrix. The constraints, Eqs. (16a, b), guarantee the Improved stiffness 
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matrix to be real symmetric. A unique set of K will result if k satisfies 
Eqs. (15), (16a, b) and also minimizes the norm 

e = ||W K (K - K a ) W K || - ||k - k A || (17) 

Def i ne 


V = e + 


m n 

. 1 , .1 
i=l j=l 


“U (+Tk - f) 1J + “*j ( * T|C • f, *j 


n n 


+ k 1 6 u <k - k *>u + 6 ij (k * - k >ij * y u < k - k >io 


+y ?j ' k - k >?j 


Note that it is no longer practical to solve the equations numerically 
for the Lagrange multipliers 3 and y because their order is n. A closed-form 
representation for the improved stiffness matrix is necessary. Differenti- 
ating v with respect to k and equating the result to zero yield 

(k - k A )* + fa + 3 - 6* + y - y T = 0 (19) 

The reason for defining Lagrangian function in the form of Eq. (18) is now 

★ 

clear since if f is differentiated with respect to k rather than k, the same 
equation, Eq. (19), is obtained. Taking the conjugate transpose of Eq. (19) 
and adding it to Eq. (19) result in 

2 (k - k A ) + <|>a + (<J>a)* + 0 + 0* = 0 (20) 

where 0 = y-y^ and the relationship k-k A = (k-k A )* has been used in the equa- 
tion. Premultiplying Eq. (20) by 4>^, applying Eq. (15), and noting that 
is nonsingular, then 

a = - U T <1>) _1 2(f-f ft ) + <(> T (0+0*)j - U 1 "*) -1 <t> T a*<f>* (21) 
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The Lagrange multipliers a and a* can be eliminated by substituting Eq. (21) 
into Eq. (20) which, after some algebraic manipulations, results in 


k = k A + g + h 

where 

g = <f> 1 (f - f A ) + j^>(<t>T<j>) l (f - f A ) 



( 22 ) 


(22a) 


+ (f - f A )* ' $ ( | t | T<f>) 


h = - j [i - <t.(4> T <t>) _1 4> T J (© + G*) I - ♦U T *r 1 4> T * (22b) 

f A = 4>\ (22c) 

Applying the symmetric condition Eq. (16b) to Eq. (22) yields 

g + h = g T + h T (23) 

Premultiply and postmultiply Eq. (23) by <f>* and respectively, and from Eq. 
(22b), note that h T <t> = 0. Then 

<|>*h<f> = <i>* (g T - g) <f> = Q (24) 

Solving Eq. (24) for 0+0* 

WjJ 1 (0 + 0*) W" 1 = - 2 ($*Z) + Q ($*Z) + j* (25) 

+ 9 

in which ($*Z) is the Moore-Penrose inverse of **Z. 


Therefore, the improved stiffness matrix becomes 

K = K A + (Z*$) + Z* * Q(Z*$) + Z* - WjJ 2 F - (W‘ 2 F)* + \ WjJ 2 PW“ 2 (26) 

where 

G = $ ($ T WjJ 2 $) _1 (27a) 

H = n 2 $ T M + fi$ T C + *\ (27b) 
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F = G H 


(27c) 

i — Z — i 
CD 
li 

Q_ 

Jh* + ($ T H*)*J g* 

(27d) 

Q = 4>* ■ 

WjJ 2 (F - F) + [w" 2 (F - F)j*+ j U ~ 2 (P T - P)W ' 2 

„ j 

- 4 (27e) 

Z = I - 

W' 2 G 4> T 

( 27f ) 


The matrix inversions required for the evaluation of the improved damping mat- 
-2 T -2 -1 

rix are and ($ 4>) . Neither will cause difficulty since the first one 

is a chosen weighting matrix and the other is of order m. 


Two important features are related to the selection of weighting matri- 
ces. First, the present identification is based on the minimum-norm approach, 
thus, the chosen weighting matrices reflect, in a sense, the confidence level 
in each element of the original analytical model. Second, the improved mass 
and damping matrices are used to identify the stiffness matrix. For large 
structural systems, it is quite impractical to specify confidence levels for 

all elements of system matrices. A convenient and physically reasonable 

- 1/2 - 1/2 

choice is to use M. and C. ' as the weighting matrices for mass and 

" " - 1/2 

damping identification and then M ' ' for stiffness identification. Other 
options, provided in the "Computer Implementation" (see page 19), are identity 
matrix and user supplied diagonal matrices. 


Note that the analysis up to Eq. (23) is exact. Numerical solution for 0 

from Eq. (23) would be prohibitive because the unknown matrix is of order n. 

★ 

Eq. (24) is obtained through pre- and post-multiplying h by $ and <(>. The 
★ 

rank of <|> h<t> in Eq. (24) is at most m (the rank of <(>) which is in general less 
than that of h. It is believed that the analysis is only approximate after 
Eq. (24) because of the reduction in rank due to the pre and post multiplica- 
tions. The inexact orthogonality check in the example are explained by this 
consideration. However, it appears that a similar procedure is exact in case 
of undamped system identification^. Further research in this area is 

required. 
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COMPUTER IMPLEMENTATION 
General Description 


Fig. 1 summarizes the computer program designed for improving mass, 
damping and stiffness matrices described in this report. A separate program, 
AMIR, is used to reorder analytical matrices on sequential files to place test 
degrees of freedom in the upper left of each matrix (see Eq. (8b)) before 
invoking the present program. 

The program performs the following functions: 

(1) Input and validate data. The data input include n, m, J l t test modal 
parameters (natural frequencies, damping ratios and mode shapes) and weighting 
matrices (optional). 

(2) Full mode computation. If m < n, Eq. (9) is solved to compute the 
full modes. A user option is available to ignore the frequency dependent 
terms (Guyan reduction) or using LDL^ (see Ref. 6) decomposition to solve Eq. 
(9) exactly. 

(3) Normalize the mode shapes by using the orthogonality relationship, 
Eq. (6a). 

(4) Formulate E and A in Eq. (14) and solve for Lagrange multipliers. 
The CPU time spent in performing this job depends on the number of measured 
modes, m, and may be the most time consuming step if m is not small. 

(5) Compute the improved mass and damping matrices by using Eq. (13) and 
check the orthogonality relationship, Eq. (6). 

(6) Compute the improved stiffness matrix by using Eqs. (26) and (27) 
and check the orthogonality relationship, Eq. (7). 
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Figure 1. Program AMIMP 










(7) The improved mass, damping and stiffness matrices are listed and 
placed on sequential files and statistical data is computed and listed. This 
data consists of: RMS of original matrix, RMS of changes, the ratio of the 
preceding, the mean absolute ratio of the diagonal changes to the original 
diagonals, the RMS of the changes divided by the corresponding diagonal ele- 
ments, i.e., the square root of the mean of 


I 


i > J 


(M 




m a } 

M jj 


1/2 


In addition, the 50 largest changes are printed. 


Fourteen sequential files are used in the program, one for the input 
data, 3 for analytical matrices (M^, C^, K^), 3 for improved matrices (M, C, 
K) and the rest are working files. At the completion of the implementation, 
the analytical matrices will be preserved and three files will contain the 
improved matrices. 


Numerical Schemes 

(1) The LDL^ decomposition used to compute full modes is the same as 
that described in Ref. 6. This algorithm performs lower-diagonal -lower trans- 
pose decomposition of a symmetric matrix and is designed for large matrices 
which do not fit in core. For details, see Ref. 6. 

(2) E in Eq. (14b) is a full matrix containing m(m+l) x m(m+l) elements. 

Let E(p,q) be the (p,q)th element of matrix E and be the (i, k, l, 

lS 1 KX'J 

j)th element of E^ s (r=l,2; s=l,2) as shown in the Appendix. Then the coeffi- 
cient matrix E is formulated as follows 

E <P.q) - (E rs ) iktj 

where 

p = (i ( i - 1 )/2 ) + j for r=l 

= (i (i-l)/2) + j + m (m + l)/2 for r=2 
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for s=l 
for s=2 


q = (k ( k- 1 )/2) + £ 

= (k (k-l)/2) + £ + m (m + l)/2 


i ,k,£ = 1 to m 


The vector A in Eq. (14b) contains m (m+1) elements. If A(p) denotes the pth 
element of A and (a r )^, (Aj).^ the (i,j)th element of A R , Aj respectively, 
then 

Mp) - (i R )ij or < A i>ij 

where 

p = (i ( i -1 )/2) + j for A(p) = (a r )... 

= (i (i-l)/2) + j + m (m + l)/2 for A(p) = (Aj)^ 

An LU decomposition algorithm is then applied to solve A from Eq. (14b). 

(3) The following algorithm (Ref. 11) is used to compute the Moore- 
Penrose inverse (A + ) of an arbitrary rectangular matrix A (rxs), r £ s. 


1. Compute C=A*A. Note that the size of C(sxs) is no larger than 
that of A(rxs). 

2. Use LU decomposition to find nonsingular matrices P and Q so 
that 


PC 2 3 4 Q = I 


o 



2 

where I is the identity matrix of the rank of C . 

3. Compute S = QI q P. 

4. Compute C + = (CS) (CS) C. 
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5. Compute A + = C + A*. 

Since (A T ) + = (A + )^, it is A^, rather than A, which must be 
used to compute A + for A (rxs) in steps 1 to 5, if r < s. 

User Options 

There are some user options in the program: 

(1) Input/Output 

ITYP = 0 normal input 

1 special NASTRAN input 
TRITE = 0 no matrix output 

1 diagonal element output, with ratio changes 

2 full matrices (analytical and changes) 

(2) Analytical damping matrix 
ICOP = 1 user supplied matrix 

2 proportional damping matrix (can be null matrix) 

If ICOP = 2, then C ft = COEM x M A + COEK x K., where COEM and COEK 
are user supplied coefficients. 
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(3) Full mode computation 


FDEP = 0 frequency independent (Guyan reduction) 

= 1 frequency dependent 
FMOD =0 full program execution 

= 1 stop when the mode computation finishes 
(4) Weighting matrices 

IWM = 0 W M = 1 

■ 1 “m 2 ■ m a 

_2 

= 2 Wjy| = user supplied diagonal matrix 

IWC = 0 w c = I 

= 1 W C 2 = C A 

_2 

=2 Wc = user supplied diagonal matrix 

IWK = 0 W K = 1 

- 1 “k 2 ■ M 

_2 

=2 = user supplied diagonal matrix 
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A NUMERICAL EXAMPLE 


A 10 DOF lumped mass model and the corresponding analytical mass, stiff- 
ness and damping matrices are shown in Fig. 2. Eigenvalues and eigenvectors 
of the system were obtained by using Kaman's computer program, DYSCO. Fig. 3 
lists computed eigensolutions of the first two modes (m=2) with five test 
points U=5). The "measured" modes n and 4> were simulated by perturbing up to 
10% the computed n and the real part of computed 4> and multiplying the imagin- 
ary part of computed $ by 100. With 0, $, M^, C A and K^, the procedure des- 
cribed in this report was then applied to improve the analytical model. The 
results are shown in Figs. 4 and 5. In Fig. 4, the diagonal elements of M^ 
and M and the percentage change of these diagonal elements are listed on the 
first three columns. Columns 4 and 5 are maximum absolute changes of the ele- 
ments beyond the diagonal element in each row and the percentage change with 
respect to the corresponding diagonal element in M^. In this example, both 
the maximum percentage change of column 3 (28%) and column 5 (33%) appear in 
the 3rd row of the mass matrix. For the damping matrix, the maximum percen- 
tage change of the diagonal elements is 8.5%. Three columns (9 to 11) which 
list the sub-diagonal element in each row of C A and C and the percentage 
change with respect to corresponding C^ element are also included. Columns 12 
and 13 give the maximum absolute change of element beyond the diagonal and 
sub-diagonal ones in each row. Similarly, percentage changes of the diagonal, 
sub-diagonal elements and the maximum change in each row are given in Fig. 5 
for stiffness matrix. 

Given in Fig. 6 is a partial validation of improved model by checking the 
orthogonality relationships, Eqs. (6a) and (7a). As shown in the table, the 

improved mass and damping matrices excellently satisfy Eq. (6a) with numerical 

-8 

errors of order 10" . The improved mass and stiffness matrices, however, do 
not satisfy Eq. (7a) as well. The resultant matrix is very close to the mea- 
sured fi except that the second damping ratio here reads -0.01414 rather than 

_3 

0.00916. The off-diagonal terms are in error by about -10 . In this regard, 

improvement of stiffnes matrix is only partially successful. 
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An eigenanalysis has been performed on the improved model. For compari- 
son, Fig. 7 lists the first two modes of the improved model, as well as the 
measured modes duplicated from Fig. 3. It is noted that the imaginary part of 
and real part of $ of the improved model are almost identical to their mea- 
sured correspondence. However, the real part of n and imaginary part of $ are 
in error. Again, this is attributable to an imperfect improvement of the 
stiffness matrix. 


?1 



30 -10 
-10 15 -5 

-5 15 -10 

-10 15 -5 

-5 15 -10 

-10 15 -5 

-5 15 -10 

-10 15 -5 

-5 15 -10 

-10 30 

' 0.6 - 0.2 

-0.2 0.3 -0.1 

-0.1 0.3 -0.2 

-0.2 0.3 -0.1 

-0.1 0.3 -0.2 

C A = -0.2 0.3 -0.1 

-0.1 0.3 -0.2 

-0.2 0.3 -0.1 

-0.1 0.3 -0.2 

- 0.2 0.6 



Figure 2. Analytical model. 




COMPUTED MODES 


MEASURED MODES 


'-0.002486 + 0.49861 


'-0.002586 + 0.5286i 


Q = 


-0.009555 + 0 . 9774i_ 


-0.009155 + 0.9374i_ 


0. 1119-1. 46E-6i 

-0. 2159-1. 6E-7i“ 


' 0. 1189-1. 46E-41 

-0. 2259-1. 6E-51' 

0.3301-3.99E-61 

-0. 6064-4. 2E-7i 


0. 3401-3. 99E-41 

-0. 6364-4. 2E-5i 

0.7008-6. 54E-6i 

-0. 9239-4. 7E-7i 

$ = 

0. 6808-6. 54E-4i 

-0. 9539-4. 7E-5i 

0.8513-6. 34E-6i 

-0. 9062-3. 9E-7i 


0. 8913-6. 34E-4i 

-0. 8762-3. 9E-5i 

0. 9830-2. 35E-6i 

-0. 1779+9. 0E-8i 


0. 9530-2. 35E-41 

-0. 1879+9. OE-61 


Figure 3. Computed and measured (perturbed) modes. 
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Figure 4. Improved mass and damping matrices. 
















































Element 

3> T M + 0 

<i> T M $ + $ T C $ 

0* T M - 

* T K 4 

Real 

Imaginary 

Real 

Imaginary 

1,1 

1.0 (1)* 

-1.3 x 10" 8 (0) 

-2.538 (-2.537) x 10' 3 

5.286 (5.286) x 10" 1 

1,2 

0 

cn 

1 

o 

H 

X 

CM 

5.8 x 10‘ 10 (0) 

3.263 x 10' 3 (0) 

8.24 x 10‘ 6 (0) 

2,1 

1.7 x 10' 8 (0) 

-1.1 x 10' 9 (0) 

3.263 x 10" 3 (0) 

8.24 x 10" 6 (0) 

2,2 

1.0 (1) 

-3.6 x 10' 9 (0) 

-1.414 (0.9155) x 10" 2 

9.374 (9.374) x 10" 1 


*Exact values are given in parentheses. 


Figure 6. Orthogonality check for the improved model 



























MEASURED MODES 


FIRST TWO MODES OF IMPROVED MODEL 



■-0. 002586 + 0.5286i 


■-0. 002537 + 0.52861 

n = 


fi = 



-0.009155 + 0. 9374i_ 


-0.01415 + 0.9374i. 


0. 1189-1. 46E-4i 

-0. 2259-1. 6E-5i' 


‘ 0. 1189-1. 29E-31 

-0 . 2258— . 791E-2i" 

0. 3401-3. 99E-41 

-0.6364-4.2E-51 


0. 3400-3. 49E-3i 

-0.6363-2.28E-2i 

0. 6808-6. 54E-41 

-0. 9539-4. 7E-5i 

$ = 

0. 6807-5. 58E-31 

-0. 9539-3. 81E-21 

0. 8913-6. 34E-4i 

-0. 8762-3. 9E-5i 


0. 8913-6. 59E-3i 

-0. 8764-4. 57E-21 

0. 9530-2. 35E-4i 

-0. 1879+9. OE-61. 


. 0. 9530-3. 81E-31 

-0. 1887-2. 30E-2i_ 


Figure 7. Measured inodes vs modes of improved model. 



CONCLUDING REMARKS 


A method is developed to find the smallest changes in the analytical 
model of a nonproportionally damped structure so that the improved model is 
compatible with the test modal parameters. The only assumptions about the 
system damping are that it is viscous and can be represented by a real symmet- 
ric matrix. The reported procedure aims at applications to large structural 
systems and has used methods which are numerically feasible. 

The following conclusions can be drawn from previous discussion: 

1) Given an incomplete set of modal parameters and an analytical model 
of a structure, application of this procedure will lead to an improved model. 
Changes between two models serve as a useful guideline for the validity of the 
original model being a good representation of the structure. The engineering 
judgement required using this approach is to establish acceptable levels of 
change. This procedure is, therefore, efficient and objective. 

(2) A numerical experiment using arbitrary perturbations in modal param- 
eters found no numerical sensitivities. 

(3) The improved stiffness matrix obtained from the present procedure, 
although better than the original one, is considered not completely satisfac- 
tory since M and C satisfy the orthogonality relationship, Eq. (6a), but M and 
K do not exactly fit the orthogonality relationship, Eq. (7a). Further 
research is required to overcome this deficiency. 

(4) The improved matrices are filled using the present procedure and may 
imply non-existent coupling. This effect may be insignificant if the resul- 
ting coupling terms are small. Additional constraints are required and fur- 
ther research is necessary to keep this problem of a solvable size. 
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APPENDIX 


Explicit Function Forms of Coefficients in Eq. (14) 

^ll^ikJtj = ^ll^iktj + ^ll^itkj 

* E 12*ik£j = (e 12^ikJlj + (1 ' 6 k)i^ (e 12^i Jtkj 

^ E 21 ^iktj = ^ e 21 ^iktj + ^ l ~ 6 ki) ^ e 2 pitkj 

^ E 22^iktj = ^ e 22^ik£j + ^ e 22^i fkj 

(e ll^iktj = 

‘ ^ M RR^ik ^ M UU^tj " ^ M RU ^ i k ( M Rlp£j + ^ M RS^ik ( M VlPtj + ^ M RV^i k ^ M SU^j 

‘ ^UR^ik " ^ M UU^i k ^RR^tj + ^US^ik ^VR^j + ^UV^Ik ( M SR^j 

+ ^SR^ik ^UV^Jlj + ^ M SU ^ i k ( M RV^j " ^SS^ik ^ M VV^j " ^SV^ik ( M SV^j 

+ (M VR ) ik ( M us^tj + ^ M VU^ik ^ M RS^tj ' ^Vik ( M VS^j “ ^Wik ^ M SS^j 

" ^ C RR^ik ^RR^Jtj + ^ C RS^ik ^ C SR^j + ^ C SR^ik ( C RS^j “ ^SS^ik ^ C SS^j 

^ e 12^iktj = 

" ( M RR^ik ^VllUj " ^Vik ^ M RU^tj " ^RS^ik ^ M UU^j " ( M RlPik ^ M SlPtj 

" ^ M UR^ik ^ M VR^tj " ^ M UV ^ i k ( M RR^j ' ^US^ik ( M UR^j ' ^ M UU^ik ^ M SR^j 

+ ^SR^ik ^ M VV^j + ^ M SV^ik ^ M RV^j + ^SS^ik ^ M UV^j + ^ M SlPik ^ M SV^j 

+ (M VR ) ik ( M vs^£j + ^ M VV ^ i k ( M RS^j + ^ M VS^i k ^US^j + (*Vik ^ M SS^j 

" ^ C RR^ik ^ C SR^tj “ ^ C RS^i k ^RR^j + ^SR^ik ^ C SS^j + ^SS^ik 
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( e 21^ik£j = 

^RR^ik ^ M UV^j + ( M RlPik ^ M RV^j " ^\shk ( M w\j " (*Wlk ^ M SV^j 
+ ( M vR^ik ( M UR^j + ^VlPik ^ M RR^j " ^VS^ik ^ M VR^j " ^VV^Ik ^ M SR^j 

+ ^SR^ik ^ M UU^j + ( M SlPik ^ M RU^j " ^ M SS^ i k ^ M VU^j " ^SV^ik ^ M SU^j 

+ ^ M UR^ik ( M US^j + ^UlPi k (M RS^j " (M US^ik ^ M VS^j ' (N Wik (M SS^j 

+ ^ C RR^ik ^ C RS^j " ^ C RS ^ i k ^ C SS^j + ^SR^ik " (C SS^ik (C SR^j 

( e 22^ik£j = 

^RR^ik ^ M VV^j + ^ M RV ^ i k ^ M RV^j + ^ M RS^ik ^ M UV^j + ( M RlPik ^ M SV^j 
+ ( M vR^ik ^ M VR^j + ^Vik (f W*j + (M VS^ik (M UR^j + (f Wik (M SR^j 

+ ^ M SR ^ i k ( M VU^j + ^SV^ik * M RlPtj + ^SS^ik ^ M UU^j + ^ M SU^ik ^ M SU^j 

+ ^UR^ik ^ M VS^j + ^UV^ik ^ M RS^j + ^US^ik ^US^j + (*Wik ^ M SS^j 

+ ^ C RR^ik ^ C SS^j + ^Wik ^ C RS^j + ^SR^ik ^ C SR^j + ^SS^ik ^ C RP^j 

where 


6 


k* 


= Dirac delta function 


and for P, Q denoting R, S, U, or V, 


m pq = P V 0 

Cpq = P T W- 2 Q 

Note that M p g = mJ r and C p g = cj 
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